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We cross correlate the gravitational lensing map extracted from cosmic microwave background 
measurements by the Wilkinson Microwave Anisotropy Probe (WMAP) with the radio galaxy dis- 
CNl 1 tribution from the NRAO VLA Sky Survey (NVSS) by using a quadratic estimator technique. We 

use the full covariance matrix to filter the data, and calculate the cross-power spectra for the lensing- 
■ galaxy correlation. We explore the impact of changing the values of cosmological parameters on 

s ! ■ the lensing reconstruction, and obtain statistical detection significances at > 3o\ The results of all 

^ ' cross correlations pass the curl null test as well as a complementary diagnostic test using the NVSS 

. data in equatorial coordinates. We forecast the potential for Planck and NVSS to constrain the 

' lensing-galaxy cross correlation as well as the galaxy bias. The lensing-galaxy cross-power spectra 

are found to be Gaussian distributed. 

CN i PACS numbers: 98.70.Vc, 98.62.Sb, 98.80.Es 
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I. INTRODUCTION 

The cosmic microwave background (CMB) tempera- 
ture anisotropy contains a wealth of cosmological infor- 
mation and has played a pivotal role in our understanding 
of the Universe. Besides the primordial fluctuations, vari- 
ous secondary anisotropies, e.g. gravitational lensing, the 
thermal Sunyaev-Zcl'dovich effect, the kinetic Sunyaev- 
Zcl'dovich effect, as well as the integrated Sachs- Wolfe 
effect, are playing an increasingly important role in con- 
straining cosmological constituents and dynamics. 

Among the secondary effects imprinted on the CMB 
gravitational lensing is of great importance. The pro- 
jected gravitational lensing potential is a line-of-sight 
probe which contains information about the geometric 
distance traversed by CMB photons and time-dependent 
gravitational potentials. As such it is very sensitive to 
late universe parameters, such as the sum of neutrino 
masses, the dark energy equation of state and spatial 
curvature. Since the projected gravitational lensing po- 
tential contains both geometric and structure growth 
information, it effectively breaks the angular diameter 
distance degeneracy [1]. Gravitational lensing measure- 
ments can also be used to de-lens the -B-mode polar- 
ization of the CMB [2], enabling us to learn about pri- 
mordial gravitational waves [3] and the energy scale of 
inflation. 

Tentative CMB weak lensing searches have been done 
with WMAP-7 year data sets [4, 5] using non-Gaussian 
statistics. However, WMAP-7 alone cannot detect weak 
lensing of the CMB because WMAP temperature maps 
have insufficient sensitivity [5]. Recently, the Atacama 
Cosmology Telescope [6] and South Pole Telescope 
(SPT) [7] have performed the first internal lensing re- 
construction detections using non-Gaussianity. In addi- 
tion, Atacama Cosmology Telescope and SPT also mea- 
sured the gravitational lensing signal from the smoothing 



effects of the acoustic peaks on the CMB temperature 
power spectrum [8-10]. As the experimental sensitivity 
improves, internal measurements, either from the power 
spectrum or the trispectrum, will become more precise 
in the near future. 

The correlation between lensing and large scale struc- 
ture arises from large scale structure, which deflects CMB 
photons in the late universe. The signal-to-noise ratio 
of lensing measurements can be enhanced if the CMB 
maps are cross correlated with highly sensitive large scale 
structure tracers, such as luminous red galaxies (LRGs) 
(which cover the redshift range 0.2 < z < 0.7), quasars 
(which covers the redshift range z < 2.7) from the Sloan 
Digital Sky Survey (SDSS), or the NVSS of radio galax- 
ies which has a higher mean redshift (z ~ 1) than the 
LRGs and quasars. Hirata et al. [11] used the cross cor- 
relation between WMAP-1 and LRGs and quasars from 
SDSS imaging and found no statistically significant sig- 
nal. Then Smith et al. [12] used the cross correlation be- 
tween WMAP-3 and NVSS, and found a 3.4er signal, in- 
cluding systematics. Using a slightly less optimal estima- 
tor than Ref. [12], Hirata et al. [13] obtained results con- 
sistent with, though at slightly lower significance than, 
Ref. [12] for WMAP-3 with LRGs (0.95cr), WMAP-3 with 
quasars (1.64cr), and WMAP-3 with NVSS (2.13er) re- 
spectively. Recently, SPT found a greater than 4a cross 
correlation between the SPT convergence field and the 
galaxy survey from the Blanco Cosmology Survey, the 
Wide-field Infrared Survey Explorer, and Spitzcr [14]. In 
this work, we use WMAP data released in years 1, 3, 5, 7, 
with NVSS to probe the lensing-galaxy correlation. We 
follow the methods developed in Smith et al. [12] and Hi- 
rata et al. [13] using all of WMAP's datasets and compare 
our results to these earlier analyses. 

The structure of the paper is as follows. We introduce 
the data sets in Sec. II. Gravitational lensing effects on 
the CMB as well as the lensing extraction technique are 
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reviewed in Sec. III. We describe the cross correlation es- 
timators in Sec. IV, and the forecast for Planck in Sec. V. 
We discuss our results in Sec. VI. 



II. DATA SETS: WMAP AND NVSS 

The CMB data we use are from WMAP's Q-, V-, and 
W-band raw differencing assemblies (DAs). All of these 
DAs are masked by the KpO mask (Fig. 1) to remove 
bright sources and the galactic plane leaving 77% of the 
sky. 

The input for the galaxy distribution is the NVSS of 
radio galaxies. The NVSS [15] team provides the software 
"NVSSlist" to convert its raw catalog to a deconvolved 
one which is corrected for known biases and systematic 
errors. We use the deconvolved catalog to extract the 
galaxy count map. We use this software here, without 
specifying either a minimum or maximum flux cut. The 
NVSS map is pixelized with a HEALPix pixelization 
scheme with A^ide = 256. We remove the galactic plane 
(\b\ < 10°) and the part of the sky unobserved by the 
survey (6 < —36.87°). We also carefully remove bright 
sources with flux > Uy and mask out a disk of radius 1° 
around them, forming the NVSS mask shown in Fig. 1. 
The resulting galaxy count map has 1, 224, 990 sources, a 
sky fraction / s k y = 0.573, a mean number of sources per 
pixel n = 2.72, and a surface density of 170,249 galaxies 
per steradian. This agrees well with previous studies [16]. 



III. GRAVITATIONAL LENSING OF THE CMB 

The effect of lensing on the CMB's primordial temper- 
ature T in direction n can be represented by 



T(n) = T(n + d(n)), 



(1) 



where T is the lensed temperature and the deflection an- 
gle field d(n) = V0, with cj> being the lensing potential. 
The operator V is the covariance derivative on the sphere 
with respect to the angular position n. We use Gaussian 
natural units with K = c = 1 throughout this paper. 

The two-point correlation function of the temperature 
field is [17]: 

(TlmTl'm') = Cf T 8ui8 m - m i{—l) m + 

Et- 1 )*' ( m L> —M ) ^ ^ LM ' (2) 

where the second term encodes the effects of lensing with 
the weighting factor ffjj, given by 



fiw — Q qFvli + Cv oFili'- (3) 
We use the standard spherical harmonic decomposition 

T(n)=^T im y im (n), (4) 



which defines the temperature modes T/ m . We use a sim- 
ilar notation for all other quantities defined on a sphere, 
e.g. (j)LM are the modes of the lensing potential, etc. 

Here Cj T is the unlensed temperature power spec- 
trum, and 



- / (2f + l)(2r + l)(2L + l) „ 
° FlLl ' = V 4^ X 

I[L(L + l) + Z'(Z' + l)-Z(Z + l)](j o o) (5) 

are proportional to the Wigner 3j symbols. Equation (2) 
provides a way to extract (j)LM from the TT correlations. 

In the late universe, the Poisson equation relates the 
lensing potential $(k) to the density contrast <5(k), 



(0) 



where Q m is the matter fraction, a is the scale factor and 
H is the Hubble constant. Using the definition D(x) = 
— 2(l/x — l/x*)i the projected lensing potential can be 
expressed as an integral along the line-of-sight, 



(n) = 



d x <f>(xn)D{ X ), 



(7) 



and it is integrated from to the comoving distance at 
the last scattering surface x*- Here x( z ) ' 1S the comoving 
distance at redshift z. The galaxy overdensity is also 
given by a line-of-sight integration as 



3(n) 



Jd X b g N( X )5( X n) 
Jd x Af( X ) 



(8) 



To better understand the projected galaxy overdensity, 
Af(x) ~ dN/dx is the comoving distance distribution 
of the galaxies. For NVSS galaxies, there is a lack of 
accurate photometric rcdshifts so approximations to the 
redshift distributions are made [18-20]. Following [12], 
we use a Gaussian distribution 



dN 

dz 



oc e 



(9) 



where a = 0.8 for z < Zq(= 1.1), and a = 0.3 for z > zq, 
and the comoving distribution dN/dx is easily derived 
fromEq. (9). 

The parameter b g (z) is the redshift-dependent galaxy 
bias. To keep the model as simple as possible we treat 
the galaxy bias as a constant which can be determined 
from a fit to the data shown in Fig. 2. This is dif- 
ferent from [13, 19] which used the cross correlation 
of NVSS with the SDSS and with sources from the 2- 
Micron All Sky Survey to determine the galaxy bias. The 
galaxy bias is of great importance because the cross cor- 
relation can be directly translated into primordial non- 
Gaussianity [21] and may enable general relativity to be 
tested on cosmological scales [22] . 

Equations (7) and (8) give the general definitions for 
the lensing potential and galaxy overdensity which are 



FIG. 1. WMAP KpO mask (left) with / sky = 0.77 and NVSS mask (right) with / sky = 0.573. 



both Gaussian fields, characterized by their variances, 
(i.e. the power spectra) Cf^ and C 99 . From Eq. (7) and 
(8), one sees that the lensing-galaxy cross correlation is 
built on the relation described by Eq. (6) . Using the Lim- 
ber approximation k ~ l/x we calculate the theoretical 
galaxy auto-power spectrum and cross-power spectrum 
in Eqs. (10) and (11). 

The galaxy-galaxy power spectrum ( (gi m 9* m ) ) is 



Cf 9 



1 



9 X 2 



l~] 



and it will be used later for determining the galaxy bias 
and for simulating galaxy maps. This power spectrum 
shows the galaxy clustering strength on different an- 
gular scales. We calculate the power spectrum of the 
NVSS overdensity map using two independent methods: 
a pseudo-C; method and a spherical harmonic estima- 
tion, as described in Refs. [18, 23] respectively. We find 
that both methods agree very well (Fig. 2). As a final 
check, we have computed the NVSS galaxy power spec- 
trum using the NVSS galaxy map in both equatorial and 
galactic coordinates and find a negligible difference, as 
expected, since the galaxy clustering is an intrinsic prop- 
erty of the Universe and does not depend on the choice 
of coordinate system. The galaxy power spectrum we 
use, obtained using the pixclized map in equatorial coor- 
dinates, is plotted in Fig. 2. 

The lensing-galaxy cross-power spectrum is 



Cf 



1 



2 Jd X M( X ) 
d X b g (l + z)D( X )^MPs(- 



k 2 x 2 



(11) 



which shows the mutual influence between the gravita- 
tional potential and the galaxy clustering in the late uni- 
verse on different angular scales. Ps(k) is the matter 
power spectrum defined using the same convention as 
Rcf . [24] . The primordial scalar curvature perturbations 
are evaluated at the pivot scale fco = 0.002 Mpc . The 
cross-power spectrum C, will be used to simulate the 
correlated galaxy maps and will also be fit to data to 
determine the detection significance. 
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FIG. 2. The NVSS galaxy auto-power spectrum. The la er- 
ror bars are from 1000 Monte Carlo simulations of the NVSS 
galaxy map with galaxy bias b g = 1. For both sets of data 
points, the red points are from the pseudo-C; method [18], 
and the green are from the spherical harmonic estimation [23]. 
The theoretical galaxy auto-power spectrum is fit to the red 
data points derived from the pseudo-C; method. Both meth- 
ods show a consistent galaxy auto-power spectrum from the 
NVSS data. The first bin of the red data points largely devi- 
ates from the theoretical curve due to systematic effects. 



IV. CROSS CORRELATION ESTIMATION 

Monte Carlo simulations are used to estimate the cross 
correlation between the CMB and the galaxy distribu- 
tion. The variances Cf T , Cf T , Cj are computed using 
CAMB [26] with the cosmological parameters listed in 
Table I. In addition to these, we derive Cf 9 and Cf 9 
from Eq. (10) and Eq. (11) with the parameters listed 
in Table I. Simulated CMB temperature modes, aim, 
are drawn from Gaussian distributions with zero means 
and variances Cf T . In this work, we will use two sets 
of cosmological parameters because we want to check the 
consistency of our results with a previous study [12] and 
also because we want to explore the impact of using the 
newest parameters from WMAP-7 on the lensing-galaxy 
cross correlations. 

We convert these a\ m to an unlensed temperature map, 
T(n), on which we do a cubic interpolation to precisely 
implement Eq. (1). This produces a lensed temperature 
map T(n) that is converted back to harmonic space to 
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TABLE I. The 6-parameter ACDM model used for the simulations of the temperature, galaxy and lensing potential. The derived 
parameter erg, based on the 6-parameter model, is shown in column eight. Using these parameters 1000 galaxy simulations 
with b g = 1 were performed to get the reconstructed galaxy biases as well as the la error bars. From column nine, we see that 
all the reconstructed galaxy biases are consistent with the input value b g — 1. Furthermore, the galaxy biases of the real data 
are calculated based on the simulations and shown in column ten. Two independent methods were used to calculate the galaxy 
auto-power spectra, as specified in the footnotes. 



Data set 


VL b h 2 




Ho 




n s 


T 


C8 


l sim 
h 3 


7 data 

°g 


WMAP-1 C 


0.0226 


0.1104 


72 


2.212 x 10" 9 


0.96 


0.117 


0.76 


0.98±0.12 a 
0.96±0.10 b 


2.00±0.12 a 
1.95 ±0.10 b 


WMAP-3 d 


0.02186 


0.1105 


70.4 


2.393 x 10" 9 


0.947 


0.073 


0.77 


0.98±0.11 a 
0.96±0.09 b 


1.97±0.11 a 
1.92 ±0.09 b 


WMAP-5° 


0.02305 


0.1182 


69.7 


2.484 x 10" 9 


0.969 


0.094 


0.85 


0.98±0.10 a 
0.96±0.08 b 


1.85±0.10 a 
1.81±0.08 b 


WMAP-7 f 


0.02258 


0.1109 


71 


2.43 x 10 -9 


0.963 


0.088 


0.80 


0.98 ±0. ll a 
0.96±0.09 b 


1.91 ±0. ll a 
1.86±0.09 b 



a Pseudo-C; method [18], 

b Spherical harmonic estimation [23] 

c WMAPl+CBI+ACBAR+2dFGRS+Lya, [25] 

d WMAP3 ALL, http://lambda.gsfc.nasa.gov/product/map/dr2/parameters.cfm 

e WMAP5+B AO+SNALL+LyoPOST, http : //lambda . gsf c . nasa . gov/pr oduct/map/dr3/par ameter s . cf m 
f WMAP7, http : //lambda. gsf c .nasa. gov/product/map/dr4/parameters . cf m 



give the lenscd modes ai m . Then each DA's beam and 
pixel window transfer function (the pixel window transfer 
function has negligible effects on the cross-power spectra) 
from WMAP are multiplied by these modes which are 
subsequently transformed into a temperature map con- 
taining the lensing signal. 

We then simulate Gaussian noise in map space where 
the pixel noise is assumed to be uncorrclated and Gaus- 
sian distributed with zero mean and pixel-independent 
variance determined from <jq / \J N b s ■ Here, both ao and 
A'obs are supplied by the WMAP team for different DAs. 
We add this noise map to the signal map and apply 
the KpO mask to get a simulated WMAP DA made in 
the same way as the real WMAP maps were produced. 
The entire procedure can be summarized by Eq. (12) in 
which ai, n is the lensed CMB mode, n(n) the white noise, 
M WMAP (n) the mask map, v the index of the DA chan- 
nel, pi the pixel window transfer function, bi the beam 
transfer function 



T^\n) = M 



WMAP 



(n) 



^2pib\^ai m Yi m (n) 

Im 



V^obs(n) 



n(n) 



(12) 



To maximize the signal-to-noise ratio, we compute a 
single harmonic mode a/ m from eight Q, V, W-band DAs. 
This reduction step is expressed as [12] 



a = (S + NO^a 
= S- 1 / 2 A- 1 S 1 / 2 N- 1 a. 



(13) 



Here a is the vector of a; m s, S the signal covariancc 
matrix, N the noise covariancc matrix, and A = I + 



g l/2 N -lgl/2 i We use the second form of Eq ( 13 ) 

and filter the raw CMB modes using a multigrid- 
preconditioned-complex conjugate gradient method. The 
master equation that has to be solved is 



Ax = y, 



(14) 



where x = S^a, and y = S 1 / 2 N~ 1 a. Equation (14) is 
better for numerical computations because A is close to 
the unit matrix. Appendix A gives details of the numer- 
ical calculation of Eq. (14). We solve Eq. (14) with ai m 
for both the temperature (T; m ) and the galaxy (gi m )- 

We use the standard quadratic estimator to recon- 
struct a noisy lensing potential map (f>i m in harmonic 
space [27, 28], 

J24>lmYi m (n) = V l ( A T {n)V l B T {n)), (15) 



Im 



where 



and 



3 B» 



E 



Q TT T/ m y Jm (n). 



(16) 



(17) 



In the above equations, Vi is the gradient operator on 
a sphere and V 1 = 5 u Vj . Here gfy is the metric of a 
sphere. 

We also use Monte Carlo simulations for the NVSS 
galaxy maps. The simulated galaxy modes g\ m are drawn 
from a Gaussian distribution and transformed into a 
galaxy overdensity map 17(11) at HEALPiX resolution 
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-Vsidc = 1024. The galaxy modes must satisfy the correct 
galaxy-galaxy auto-correlation and lensing-galaxy cross 
correlation. From these two constraints the simulated 
galaxy mode must be 



a 



gim 



cf 



cf 9 



cf 



-Gh- 



(18) 



where Gi m is a complex Gaussian random variable, and 
4>im is inherited from the deflection field in Eq. (1). From 
this equation we see that the lensing-galaxy correlation 
is encoded in the first term. A NVSS map is generated 
where the galaxy number count in each pixel is drawn 
from a Poisson distribution with mean 



A(n) = n(l + g(n)). 



(19) 



The galaxy count map A(n) is used to generate a simu- 
lated galaxy overdensity map g( slm ^(n), 
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(sim) 



(n)=M NVSS (n) 



A(n) 



1 



(20) 



where M NVSS is the NVSS mask shown in Fig. 1. 
g( slm )(n) automatically contains the shot-noise with the 
variance Nf 3 = 1/n for the galaxy overdensity map. We 
degrade this map to resolution N s a e = 256 i.e. the same 
as the real NVSS data. The harmonic mode fffc , which 
contains the shot-noise, is obtained from g( slm >(n) and is 
further filtered using the same procedure as in Eq. (13), 



ra i tvt\— 1 (sim) 
gi m = (S + N) g\ m '. 



(21) 



Here S represents the primordial galaxy covariance and 
N is the shot-noise covariance. 

We show the noisy reconstruction of the potential maps 
and the filtered galaxy map in Fig. 3, using the measured 
WMAP and NVSS data. 

The lensing-galaxy cross-power spectrum is the observ- 
able which will be compared with the counterpart from 
data. The estimator of the lensing-galaxy cross correla- 
tion is expressed as 



CI 



1 



E 

leb 

-Km<l 



{<j>lm - {4>lm))* fan 



(22) 



where JFb is the normalization factor at band b. It is 
shown in Rcf. [29] that the normalization factor can be 
calculated by either a direct or a fast method for the 
full-sky coverage and that the fast method converges in 
O(10 2 ) simulations. When there is a sky-cut these meth- 
ods account for the sky-cut effect very well and a constant 
/sky is often used. The factor / s k y is actually a function 
of I [30] so a simple constant approximation may poten- 
tially bias the cross-spectra reconstruction. Therefore an 
end-to-end simulation [12] is the best way to get the ex- 
act normalization accounting for the sky-cut and that is 
done here. 



As a systematic check we note that the lensing signal 
consists of a gradient and a curl component [31]. The 
curl component estimator tpi m is defined by 



I in 



i>lmYi m {n) =J2^i( oA T (n)V, B T (n)) (23) 



and should vanish because lensing does not generate vor- 
ticity. Similar to Eq. (22), the curl-galaxy cross correla- 
tion diagnostic is calculated by 

Cf 9 = ^ J2 (i>lm-(tPlm)Ygim (24) 

^ 6 leb 
-l<m<l 

which should also vanish. 

The amplitude of the cross correlation is determined 
using 



C 



^2ab C 'a 



(th) 



r<— i /^f( obs ) 

^AB^B 



l^AB A 



(th) 



^AB^B 



(th) 



(25) 



Cab is the covariance matrix for the band powers and 
A and B stand for the band power index. We find that 
the off-diagonal correlations of Cab are negligible, and 
the covariance matrix elements can be simply replaced 
by the band power variance a\, i.e. Cab = <*\5ab- 

We have described the procedures used to perform 
analysis on simulated or measured data. Now we sum- 
marize the analysis of the real WMAP and NVSS data. 
We fit the theoretical galaxy auto-power spectrum to the 
NVSS data in Fig. 2 and determine the galaxy biases 
(Table I) using two methods. The error bars are de- 
termined from 1000 simulated galaxy maps with galaxy 
bias b g = 1. Then we choose two sets of parameters (la- 
beled "Set-3" and "Set- 7") in Table II to do the cross 
correlation calculations in this work. In order to com- 
pare our results with those from the previous studies [12], 
we use the parameters they used Set-3 from WMAP-3's 
cosmological parameters ( row "WMAP-3" in Table I) 
combined with the corresponding galaxy bias of Smith 
et al. [12]. Based on the newest cosmological parameters 
from WMAP-7 ( row "WMAP-7" in Table I) we con- 
struct a new parameter set Set-7 with the corresponding 
galaxy bias shown in Table I. For each of the parameter 
sets we calculate four lensing-galaxy cross correlations 
from WMAP-1 to WMAP-7. 

We carefully treat the known systematics of NVSS, 
i.e. the gradient stripes which are generated by the 
declination-dependence of the galaxy overdensity field 
due to the low- flux calibration issue [15]. We first make 
a gradient stripe map only using m = modes and 
then subtract it from the galaxy map. We calculate 
the lensing-galaxy cross correlations for two cases: with- 
out the gradient stripes removed and with the gradient 
stripes removed. We find that this systematic effect does 
not affect the lensing-galaxy cross correlations as seen 
from column C a and column C b in Table III and Ta- 
ble IV. The statistical results are those with the gra- 
dient stripes removed which are shown in Figs. 4 and 5 



FIG. 3. The noisy reconstruction of the lensing potential map (Eq.(15) from WMAP-7) band-pass filtered from 20 < I < 40 
(left). The analogous map from NVSS galaxy data [Eq.(21)] within the band 20 < I < 40 (right). 



TABLE II. The two sets of cosmological parameters used in this work: we choose two sets of parameters (labeled "Set-3" and 
"Set- 7") to do the cross correlation calculations in this work. In order to compare our results with those from the previous 
studies [12], we use the parameters they used Set-3 from WMAP-3's cosmological parameters ( row "WMAP-3" in Table I) 
combined with the corresponding galaxy bias of Smith et al. [12] . Based on the newest cosmological parameters from WMAP-7 
( row "WMAP-7" in Table I) we construct a new parameter set, Set-7 with the corresponding galaxy bias shown in Table I. 



Data set 


Q. b h 2 


£2cdm/i 2 


Ho 


A s 


n s 


T 


0"8 


b g 


Set-3 


0.02186 


0.1105 


70.4 


2.393 x 10~ 9 


0.947 


0.073 


0.77 


1.70 


Set-7 


0.02258 


0.1109 


71 


2.43 x 10 -9 


0.963 


0.088 


0.80 


1.91 



for the KpO+NVSS mask combination. From the two fig- 
ures, we find that the lensing-galaxy cross-power spectra 
are consistent with the theoretical predictions and the 
uncertainty of the cross-power spectrum is decreasing as 
the year of WMAP increases. All the error bars are cal- 
culated from 1000 Monte Carlo simulations, which we 
confirmed to be sufficient for convergences. As a com- 
plementary diagnostic test, we keep the NVSS galaxy 
overdensity map in equatorial coordinates and calculate 
the cross-power spectra and all the amplitudes arc shown 
in column C c in Table III and Table IV. As can be seen, 
they are negligible. All the cross correlation amplitudes 
are summarized in Table III and Table IV. From the re- 
sults of WMAP-3 x NVSS in Table III: for the statistical 
results, we get lensing detection significance level of 3.60er 
and [12] got 4.02cr. Both analyses agree quite well. We 
find the cross-power spectra from WMAP-5 and WMAP- 
7 clearly and firmly show the lensing-galaxy correlation 
at > 3<7 level for both cases. All the results are within 
the optimal bounds shown in Table V. 

Assuming there is no cosmological parity violation the 
curl-galaxy correlation should be consistent with zero. 
We show the results of the curl null tests in Figs. 6 and 
7. As expected, all the correlations are consistent with 
zero. The amplitude as well as the significance are given 
in Table VI. 

We pixelized the NVSS catalog with different 
HEALPix resolutions (e.g. iV sidc = 512, 1024) to probe 
the possible pixel artifacts that could afflict the cross- 
power spectra and because we want to examine the im- 
pact of possible long range spatial correlations. However, 
we find that different NVSS pixclization resolutions do 



not affect the cross correlation. 

We also use the diagonal elements of the covariance 
matrix to do the analysis to check the consistency with 
previous studies. In this case, the estimator has a larger 
variance as pointed out by Smith et al. [32]. This con- 
tributes to the difference in significance levels between 
4cr [12] and 2a [13]. 



V. FORECAST FOR FUTURE EXPERIMENTS 

The revealed cross correlation between WMAP and 
NVSS hints that the detection significance would be fur- 
ther enhanced if the precision of the CMB data were 
improved. The upcoming Planck data will improve upon 
WMAP, so we expect that the cross correlation between 
Planck and NVSS will be more significant. To predict 
the optimal bound on the detection signal-to-noise ratio 
for lensing-galaxy cross correlation we first calculate the 
equivalent noise from the following equation 



Nf 



» N 99 



1/2 



(26) 



where Nf is the lensing reconstruction noise [17] and 
Nf 9 is the galaxy shot-noise. The efficient algorithm for 
calculating 2Vj is given in Refs. [2, 5]. This reconstruc- 
tion noise can be minimized by combining different CMB 
channels and the minimum noise is 



N, 



1 



(27) 
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TABLE III. Measure of lensing-galaxy cross correlation C and its significance C/AC using Set-3. For five columns of this 
table: the second column shows the simulation results, the third column is the case without gradient stripes removed, the 
fourth column is the case with gradient stripes removed (this column shows the statistical results of the lensing-galaxy cross 
correlations). The fifth column is the case by setting the NVSS map in equatorial coordinates as a complementary diagnostic 
test. 



Data set 


C sim C/AC 


C a C/AC 


C b C/AC 


C c 


WMAP-lxNVSS 


1.00 ±0.47 2.13a 


1.25 ±0.47 2.66cr 


1.24 ± 0.47 2.64a 


0.26 


WMAP-3xNVSS 


1.00 ±0.35 2.86cr 


1.20 ±0.35 3.43(7 


1.26 ± 0.35 3.60cr 


0.17 


WMAP-5 x NVSS 


1.00 ±0.31 3.23cr 


1.24 ±0.31 4.00cr 


1.27 ±0.31 4.10(7 


0.23 


WMAP-7xNVSS 


1.00 ±0.30 3.33a 


1.14 ±0.30 3.80(7 


1.16 ± 0.30 3.87a 


0.15 



a Without gradient stripes removed. 

b With gradient stripes removed. 

c Galaxy map in equatorial coordinate. 



TABLE IV. Measure of lensing-galaxy cross correlation C and its significance C/AC using Set-7. The format of this table is 
the same as Table III. 



Data set 


C sim C/AC 


C a C/AC 


C h C/AC 


C c 


WMAP-lxNVSS 


1.00 ±0.41 2.44cr 


1.01 ±0.41 2.46cr 


1.00 ±0.41 2.44a 


0.20 


WMAP-3 x NVSS 


1.00 ±0.31 3.23cr 


0.96 ±0.31 3.10cr 


1.01 ±0.31 3.26a 


0.13 


WMAP-5 x NVSS 


1.00 ±0.28 3.57a 


0.98 ±0.28 3.50cr 


1.01 ± 0.28 3.61cr 


0.18 


WMAP-7xNVSS 


1.00 ±0.26 3.85cr 


0.92 ±0.26 3.54cr 


0.93 ± 0.26 3.58a 


0.11 



a Without gradient stripes removed. 

b With gradient stripes removed. 

c Galaxy map in equatorial coordinate. 



TABLE V. Fisher matrix analysis for WMAPxNVSS cross 
correlation. The la error bars are determined from Eq. (28). 
We calculate two sets of the optimal bounds for this work, 
based on two sets of parameters: Set-3 (column two); Set-7 
(column three). 



Data set 


^optimal a q j f\C. 


c optimalb C / Ac 


WMAP-lxNVSS 


1±0.46 2.17(7 


1±0.39 2.56a 


WMAP-3 x NVSS 


1±0.29 3.45a 


1±0.25 4.00a 


WMAP-5 x NVSS 


1±0.25 4.00(7 


1±0.21 4.76a 


WMAP-7xNVSS 


1±0.22 4.55a 


1±0.19 5.26a 



a WMAP-3 year cosmological parameters and b g = 1.70. 
b WMAP-7 year cosmological parameters and b g = 1.91. 



show the signal-to-noise ratio for Planck with NVSS as 
a function of Z max in Fig. 8. We find that the highest 
signal-to-noise ratio, i.e. 15a, saturates at i max = 2000. 
Since the lensing-galaxy cross-power spectrum scales as 
Cf 9 oc b g as illustrated by Eq. (11), the amplitude of 
this cross-power spectra is degenerate with the galaxy 
bias and the signal-to-noise for the cross-power spectrum 
can also serve as a prediction of the detection significance 
for the galaxy bias. Thus, we see Planck can detect b g 
with high precision which will lead to a better under- 
standing of the correlation between the baryonic matter 
distribution and the dark matter distribution. 



Both of the noise spectra effectively propagate the uncer- 
tainty ACf 9 into the cross-power spectrum Cf 9 . Specif- 
ically, we express it as 



ACf 9 = 



(28) 



The optimal bound is then determined from 



E 



a 



ACf 



1/2 



The redshift distribution Eq. (9) was used and the 
galaxy bias was set equal to b g = 1. The instrumen- 
tal properties for Planck are given in Rcfs. [33, 34]. We 



VI. CONCLUSION 

We have calculated the lensing-galaxy cross-power 
spectra using WMAP and NVSS and the full covariance 
matrix to filter the data sets. Specifically, we performed 
a thorough analysis of WMAP-1, WMAP-3, WMAP-5 
and WMAP-7 raw DAs. The cross correlations between 
WMAP-5, -7's 8 DAs (2Q-bands±2V-bands±4W-bands) 
with NVSS clearly and firmly show signals at > 3cr level. 
We took the effects of gradient stripes into account for 
the NVSS data, and determined the significance with- 
out and with gradient stripes removed. The major ef- 
fects caused by the stripes can be seen from the first bin 
of either the galaxy auto-power spectrum or the lensing- 
galaxy cross-power spectrum; the first bin decreases if the 
gradient stripes are marginalized over. However, gradient 
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TABLE VI. Results of the curl null tests for WMAPxNVSS cross correlation. The curl null tests are performed based on two 
sets of parameters: Set-3 (column two); Set-7 (column three). 



Data set 


C a C/AC 


C b C/AC 


WMAP-lxNVSS 


-0.11 ± 0.47 -0.23a 


-0.03 ±0.41 -0.07a 


WMAP-3 x NVSS 


0.00 ±0.35 O.OOcr 


0.04 ±0.31 0.13a 


WMAP-5xNVSS 


0.05 ±0.31 0.16cr 


0.07 ±0.28 0.25cr 


WMAP-7xNVSS 


-0.05 ± 0.30 0.17a 


-0.03 ±0.26 -0.12a 



a WMAP-3 year cosmological parameters and b g = 1.70. 
b WMAP-7 year cosmological parameters and b g = 1.91. 



TABLE VII. Gaussianity diagnostics for the probability distribution of {C} which is constructed from 1000 Monte Carlo 
simulations. The second column is the Kolmogorov-Smirnov test, and the critical value is 0.04 at 5% confidence level. The 
Kolmogorov-Smirnov test requires the maximum deviation be < 0.04 to validate the distribution is Gaussian. The third column 
is the skewness of {C}, and the fourth column is the kurtosis of {C}. The upper values in the cells are the results for Set-3, the 
lower values for Set-7. For a Gaussian distribution, the skewness should be and the kurtosis should be 3. As can be seen, all 
the probability distribution functions pass the Kolmogorov-Smirnov test and are consistent with being Gaussian-distributed. 



Data set 


maximum distance[<0.04] 


skewness [~0] 


kurtosis [~3] 


WMAP-lxNVSS 


0.02 


0.02 


2.89 




0.02 


-0.05 


2.77 


WMAP-3 x NVSS 


0.02 


-0.14 


2.62 




0.03 


-0.17 


2.58 


WM AP-5 x NVSS 


0.03 


-0.17 


2.53 




0.03 


-0.23 


2.43 


WMAP-7xNVSS 


0.03 


-0.21 


2.44 




0.03 


-0.21 


2.36 



stripes do not affect the lensing-galaxy correlation (com- 
pare Refs. [20, 35, 36]). We have explicitly shown all these 
results in Tables III and IV. In these two tables, column 
C a are the results without the gradient stripes removed 
and column C b arc our main results with the gradient 
stripes removed. In order to validate the lensing-galaxy 
cross correlations, we produced a NVSS galaxy map in 
equatorial coordinates directly from the NVSS catalog 
and cross correlated it with the WMAP DA which is in 
galactic coordinates, we find that all the lensing-galaxy 
cross correlation amplitudes arc negligible. 

We investigated the impact of different NVSS pixeliza- 
tion resolutions and found no effect. We compared the 
sensitivities of the estimators both with the full and 
diagonal covariance matrix and found that the former 
more effectively reduces the variance, which is mainly 
caused by the sky-cut and the inhomogeneous instru- 
mental noise. However, the former scheme involves the 
inversion of a large matrix which is computationally chal- 
lenging. 

We predicted the detection significance for the lensing- 
galaxy cross correlation or the galaxy bias for the up- 
coming Planck data with NVSS and found the detection 
significance will be improved by a factor of 5. 

The minimum variance of the estimator assumes that 
the CMB and galaxy overdensity modes are Gaussian. 
However, if the CMB contains gravitational lensing, the 
bispectrum (TTg) is not zero; it induces an additional 



variance as indicated in Eq. (B12). We analytically and 
numerically confirm that this variance is actually not no- 
ticeable for WMAP and NVSS as pointed out in Rcf. [37] . 
Furthermore, being aware of the potential non-Gaussian 
shape of the probability distribution function (PDF) [38] , 
we specifically investigate the PDF of the cross-power 
spectrum amplitude C in terms of Kolmogorov-Smirnov 
test, the skewness and the kurtosis. The diagnostic 
tests are shown in Table VII. All the PDFs pass the 
Kolmogorov-Smirnov tests. All the PDFs are consistent 
with being Gaussian-distributed (Figs. 9 and 10). 

The lensing-galaxy cross correlations effectively link 
the early universe to the late universe and the CMB is 
served as a back light casting the dark cosmic web (which 
is formed by the dark matter) throughout the major ex- 
pansion history of the universe. The gravitational lensing 
is a powerful tool to decode the information of dark mat- 
ter distribution from the CMB and the lensing-galaxy 
cross-correlations further unveil the relationship between 
baryonic matter and dark matter. 
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FIG. 4. (Set-3) The lensing-galaxy cross-power spectra for WMAPx NVSS are calculated from Eq. (22). The KpO mask is 
used to remove the contaminated regions of the WMAP data. WMAP's data are provided from two Q bands, two V bands and 
four W bands. The NVSS mask is applied to the galaxy map to remove bright sources and unobserved regions. The theoretical 
cross-power spectra are shown in blue solid lines, and they are the same for all of the four panels. The real data are shown 
in the red scattered points. The statistical amplitude for WMAP-lxNVSS is 1.24 ± 0.47, for WMAP-3xNVSS is 1.26 ± 0.35, 
for WMAP-5xNVSS is 1.27 ± 0.31, for WMAP-7xNVSS is 1.16 ± 0.30. All the error bars are determined from 1000 Monte 
Carlo simulations. We find that the lensing-galaxy cross-power spectra are consistent with the theoretical predictions and the 
uncertainty of the cross-power spectrum is decreasing as the year of WMAP increases. 
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Appendix A: MULTIGRID-PRECONDITIONED 
COMPLEX CONJUGATE GRADIENT 
INVERSION 

Given the signal covariance matrix S and the noise co- 
variance matrix N, and an array of the CMB modes a, we 



1 http://healpix.jpl.nasa.gov/ 



define another covariance matrix A = I + S 1 ' 2 N _1 S 1 ' 2 , 
and two vectors x = S 1 / 2 a, and y = S 1 / 2 N _1 a. For 
the problem Ax = y, we write down the equations for 
constructing the matrix A and the vector y, 



N. 



Iml'm 1 



x J dn *im(n)3W(n) 



M(n) 



,(A1) 



x J dn YC m {n) 

W lm = J dn Y lm( n ) 



M(n)H(n) 



M(n) 



(A2) 
(A3) 
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FIG. 5. (Set- 7) The lensing- galaxy cross-power spectra for WMAPx NVSS are calculated from Eq. (22). See Fig. 4 for 
detailed descriptions. The statistical amplitude for WMAP-1 xNVSS is 1.00 ± 0.41, for WMAP-3 x NVSS is 1.01 ± 0.31, for 
WMAP-5 x NVSS is 1.01 ± 0.28, for WMAP-7xNVSS is 0.93 ± 0.26. 



In the above equations, pi is the window transfer func- 
tion, is the specific beam transfer function corre- 
sponding to the DA of WMAP, and M(n) is the mask 
map. For WMAP, v = Q u Q 2 , V\, V 2 , Wi, W 2 , W 3 , W±, 
H(n) = T(n) and M(n) is the KpO mask. For NVSS, 
v = 1 and H(n) = g(n). Since NVSS has 45 arc-second 
FWHM resolution [15], we set = 1 as NVSS's beam 
transfer function. The correspondence between the con- 
tinuum and discrete forms of integration on the sphere 



/ 



47T 

NT 



pix 



(A4) 



where j denotes the pixel index according to the 
HEALPix pixelization scheme and iV p i x is the total num- 
ber of pixels. 

For comparison, we also use a suboptimal estimator 
which only takes the diagonal elements of the inverse 
noise matrix NZ^, , shown in Eq. (Al). 

The filtering using the covariance matrix requires us 
to solve the linear equation Ax = y. We use the precon- 



ditioned conjugate gradient iteration to solve it, and the 
initial condition is chosen to be 

x(°) = 0, 
r (0) = y) 

with the preconditioner defined as 
A" 1 = 



Aq 1 



A- 7 

here Aa is the diagonal element of the matrix A. 
The iteration procedure is [11] 

x w = x (-d + £ — ; A , t y o, 



(A5) 



(A6) 



rW = y-AxW, 
P 



>W = A r 



i-H-i) , r A r -(i-i). (A7) 



r(*- 2 )A- 1 r(''- 2 ) 

From Eq. (A7), we find that two operations A _1 r and 
Ap arc computationally demanding if we evaluate them 
directly because A and Ap are 10 6 x 10 6 matrix. 
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FIG. 6. (Set-3) The curl null tests for WMAPx NVSS are calculated from Eq. (24). The KpO mask is used to remove the 
contaminated regions of the WMAP data. WMAP data are provided from two Q bands, two V bands and four W bands. The 
NVSS mask is applied to the galaxy map to remove bright sources and unobserved regions. The theoretical lensing- galaxy 
cross-power spectra with both WMAP and NVSS in galactic coordinates are shown in blue solid lines for comparison, and 
they are the same for all of the four panels. The curl amplitude for WMAP-1 x NVSS is —0.11 ± 0.47, for WMAP-3xNVSS is 
0.00 ± 0.35, for WMAP-5xNVSS is 0.05 ± 0.31, for WMAP-7xNVSS is -0.05 ± 0.30. As can be seen, all cross-power spectra 
for the curl null test are consistent with zero (black dotted line). 



In order to achieve the necessary efficiency, we recur- 
sively precondition the matrix A on a much coarser grid. 
The preconditioner is 



A" 1 



(A8) 



express the smallest inversion problem as follows 2 



EE 



x (-1) 




{v)a j{2h + 1)(2 Z 2 + 1)(2Z 3 + 1) 

47T 

h h h \ 

mi m 2 m 3 / 



(A9) 



then the problem of preconditioning A with Ao on the 
finer grid can be iteratively solved by using Eq. (A7). 



and on the coarser grid the preconditioner is A^ . This 
multigrid strategy enables us to directly store the matrix 
Aq on the coarsest grid and we can further analytically 



2 In the following, we denote subscripts li or £; , m; by i for sim- 
plicity, so that p t . -> p t , Y t . m . -> Yi, A^m^ma ctc ' 
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FIG. 7. (Set-7) The curl null tests for WMAPx NVSS are calculated from Eq. (24). See Fig. 6 for detailed descriptions. The 
curl amplitude for WMAP-lxNVSS is -0.03±0.41, for WMAP-3xNVSS is 0.04±0.31, for WMAP-5xNVSS is 0.07±0.28, for 
WMAP-7xNVSS is —0.03 ± 0.26. As can be seen, all cross-power spectra for the curl null test are consistent with zero (black 
dotted line). 



For this work we use three levels of the grids: (1) 
N sidc = 512, Z max = 1000, (2) N sidc = 256, Z max = 400, 
(3) A^sido = 128, Z max = 200. We split the covariancc 
matrix on the third grid at Z S piit = 30 to construct the 
minimum inversion problem. 

For the coarsest grid, we explicitly calculate the inverse 
noise matrix N^ 1 [Eq. (A9)] using LAPACK. The iter- 
ation process also needs the multiplication for AA, and 
this can be computed efficiently by doing spherical har- 
monic transformations: 



AA = ^(/ + 5 1 / 2 iV- 1 5 1 / 2 ) 1 4A4 

4 

x (E^i l/) 5 4 1/2 A 4 n U (A10) 



Appendix B: NON-GAUSSIANITY 

There are several possible non-Gaussian effects gener- 
ated by using a nonzero bispectrum in the simulation. 
These can potentially bias our results. We analytically 
calculate this non-Gaussian bias in this appendix. 

We define the bispectrum by 



(aia 253 ) = 6i 23 G(123), 



(Bl) 



where 6 123 = (/i 23 C 2 TT + foaCf^Cg 9 , (see footnote 2 
for notation) and 



G(123) = 



'(2ii + l)(2i 2 + l)(2Z 3 

47T 




( h h h 

\ mi m 2 m 3 
The estimator is 

1 



C=-(C A -C B ) 



(B2) 



(B3) 
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FIG. 8. The signal-to-noise ratio for the lensing-galaxy cross 
correlation between Planck and NVSS as a function of the 
maximum multipole used in the analysis. 

where 

Oa = IY1 ^i2 3 G(123)a 1 a 2 33 (B4) 

123 

and 

Cb = \Y< & i23G(123) [C tt ] ~l 33. (B5) 

123 

We define 

5 = C~ 1 a, 

/ fc = i£> 2fe G(12fc) [G TT ]~ 2 \ 
(aia 2 )=[C TT ]- 2 1 , 



The summation index i denotes a sum over liirii. 
We define the normalization as 



F= 2 E 6i2 3 fo456G(123)G(456) 

123456 

x^u'rtrt- ( B6 ) 



The variance of the estimator G is cr 2 (C) which has 
contributions from three parts, 



a 2 (G) = a 2 (C A ) + a 2 (C B ) - 2* 2 (C A C B ). (B7) 



Now we explicitly determine the three variances. For the 
second term, we have the relation 



1'2'3' 

(B8) 

so the last two variance terms can be easily expressed as 



<T 2 {C B ) = o 2 {C A C B ) = f T [C^y 1 f. (B9) 



For the first term, it is 



6i23^456G(123)G(456)(aia 2 g3a4a556) ~ 7 £ ^i23&456G(123)G(456)(aia 2 3 3 )(a4a 5 3 6 ), 



123456 



123456 



(BIO) 



and can be expanded as 

j £ 6 i23&456G(123)G(456)(aia 2 g 3 a4a 5 g 6 ) 



123456 



2 £ & 123^456G(123)G(456) j (ai^^jaiab^ + i"-^^^^^ + (aiQ4g3)(52Q5g6) + (Q2Q5g3)(QlQ4.96) 
123456 ^ 



second term 



(aia 5 <? 3 ) (a 2 a 4 5 6 ) + (a 2 a 4 g 3 ) (a,id 5 g e ) 



3-3 



(aia 5 ) (0,20,4) (g 3 g 6 ) + (0104) (d 2 a 5 )(g 3 g 6 ) +(0^2) (0405) (g 3 g 6 ) 



normalization 



2-2-2^ 

(Bll) 
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FIG. 9. (Set-3) Probability distribution function for the lensing-galaxy cross correlation. The likelihood functions are normal- 
ized to 1. From 1000 simulations, a set of {C} is generated for each one of the subfigures, then by counting the frequency of C 
within a bin, a step-like function (red) is plotted. For comparison, Gaussian likelihood (green) is plotted using the mean and 
the variance of the set {C}. 



a 2 (C A - C B ) = J 7 + I j ^2 6i23&456G(123)G(456) (a 4 a 5 g 3 ) (aifafa) + (61&4.93) (6 2 6 5 5 6 ) + (6265.93) (6ia 4 <? 6 ) 



123456 

+ (6165(73) (626456) + (626433) (6i6 5 g 6 ) 
0{b 2 ) + 0{b A ) 



nonzero bispectrum 



(B12) 



I 

When all the pieces are put together, we find that the 0(b 4 ). We conclude that for WMAP, this contribution is 
nonvanishing bispectrum induces an extra term which is very small and we numerically verified that this is indeed 

the case. 
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